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Using a high-statistics, high-purity sample of //^-induced charged current, charged pion events 
in mineral oil (CHz), MiniBooNE reports a collection of interaction cross sections for this process. 
This includes measurements of the CC7r + cross section as a function of neutrino energy, as well as 
flux-averaged single- and double-differential cross sections of the energy and direction of both the 
final-state muon and pion. In addition, each of the single-differential cross sections are extracted as 
a function of neutrino energy to decouple the shape of the MiniBooNE energy spectrum from the 
results. In many cases, these cross sections are the first time such quantities have been measured 
on a nuclear target and in the 1 GeV energy range. 

PACS numbers: 13.15.+g, 25.30.Pt 
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Charged-current charged pion production (CC7T+) is a 
process in which a neutrino interacts with an atomic nu- 
cleus and produces a muon, a charged pion, and recoiling 
nuclear fragments. An understanding of CC7r + interac- 
tions is important for the next generation of accelerator- 
based neutrino oscillation experiments. The signal mode 
for these experiments is the charged-current quasielas- 
tic process (CCQE), and in the few GeV neutrino en- 
ergy range where such searches are typically conducted, 
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the dominant charged-current background is from CC7r + 
events. If the pion produced in a CCtt + interaction is 
lost, the final state will be identical to that of a CCQE 
event. 

Further complicating these measurements is the use 
of nuclear targets. The spectrum of nuclear resonance 
states that most often produce the pions in CC7r + events 
is modified inside the nucleus, and interactions that pro- 
duce nucleons below the nuclear Fermi momentum are 
inhibited by Pauli exclusion. After the initial interac- 
tion takes place, the final-state particles can interact 
within the nuclear medium to absorb or produce pions, 
thus modifying the observed particle composition. When 
taken in concert, these effects make difficult the extrap- 
olation of previous measurements on hydrogen and deu- 
terium to heavier nuclei. 

Many theoretical calculations exist that predict cross 
sections and kinematics for CC7r + interactions on nuclear 
targets [IHII , but currently there are not many neutrino- 
based measurements with which these models can be 
evaluated. Measurements of the CCQE to CC7r + cross 
section ratio are available from MiniBooNE [9] , K2K [l(| , 
and ANL [llj. The only absolute cross section measure- 
ments in the 1 GeV range were conducted decades ago on 
hydrogen and deuterium bubble chambers at ANL [Til — 
HU and BNL Q. The results were based on less than 
4000 events combined and differed from one another in 
normalization by ^20%. 

MiniBooNE has collected what is currently the world's 
largest sample of CC7r + interactions: a total of 48 322 
candidate v^-CGn + events with 90% purity. A recon- 
struction algorithm has been developed to distinguish 
muons and pions based on the presence of hadronic inter- 
actions with an 88% success rate. Using this kinematic 
information, we report a measurement of the CC7r + cross 
section as a function of neutrino energy, as well as the first 
single- and double-differential cross sections for the final- 
state muon and pion. These results provide the most 
complete information available on this process as mea- 
sured on a nuclear target in the 1 GeV energy range. 



A. Signal Definition: Observable CC-n 

After the initial neutrino interaction takes place, the 
resulting final-state particles must traverse the remainder 
of the nucleus before they can be detected. This results 
in additional interactions with the nuclear medium that 
can produce or absorb pions. Since these intranuclear 
processes are not experimentally accessible, an "observ- 
able CC7r + " interaction has been defined in this analysis 
as any event with a fi~ and a 7r + leaving the nucleus, 
regardless of which particles were produced in the ini- 
tial neutrino interaction. In principle, it is possible to 
use the Monte Carlo simulation to correct the observed 
event distributions back to the initial neutrino-nucleon 
interaction, however such a correction introduces a sig- 
nificant amount of dependence on the chosen final-state 



interaction model. To reduce this model dependence, the 
measurements are reported for observable CC7r + interac- 
tions. Apart from a muon and a single pion in the final 
state (no other mesons), no requirement is made on the 
number of photons, nucleons, and multinucleon states. 

In MiniBooNE, CC7T+ interactions are dominantly pro- 
duced either through an intermediate resonance state or 
by scattering off of the entire nucleus coherently. In the 
former case, the neutrino interacts with a single nucleon, 
producing a resonance state (usually a A at MiniBooNE 
energies), which then decays to a nucleon and a pion. 
The results to follow are all combined measurements of 
both incoherent and coherent processes. 



B. MiniBooNE 

The Mini-Booster Neutrino Experiment (MiniBooNE) 
was designed to search for the appearance of oscillated 
v e events from a high-purity beam. The Booster Ac- 
celerator at the Fermi National Accelerator Laboratory 
(Fermilab) provides a beam of 8 GeV kinetic energy pro- 
tons, which is directed onto a 71 cm beryllium target. 
Positively charged particles produced in the target are 
forward-focused by a cylindricaUy symmetric horn sur- 
rounding the target. Downstream of the horn is a 50 m 
drift pipe in which the particles produced in the target 
are allowed to decay to neutrinos. These decays are dom- 
inated by 7r + — >fi + with a small contribution from 
muon and kaon decay channels. The result is a neutrino 
beam composed of 93.6% 5.9% D^, and a small con- 
tribution of v e and v e . At 1 GeV where the observed 
CC7r + event distribution is peaked, the component of 
the beam is 97% of the total flux. A detailed description 
of the beamline and the neutrino flux prediction is given 
in Ref. [l|. 

The MiniBooNE detector, located 541 m downstream 
of the target, consists of a spherical tank 610.6 cm in 
radius with a 575 cm radius main volume surrounded 
by a outer veto region, which is used to detect particles 
entering and exiting the main volume. The inside surface 
of the main volume is lined with 1280 photomultiplicr 
tubes (PMTs), and an additional 240 PMTs are mounted 
inside the veto region. Charged particles are detected 
via the light they emit as they traverse the 818 tons of 
mineral oil residing in the tank. More information on 
the performance of the MiniBooNE detector, including a 
description of the optical properties of the oil, is given in 
Ref. 0. 

C. Neutrino Interaction Simulation 

The MiniBooNE Monte Carlo simulation uses the NU- 
ANCE event generator to simulate neutrino interac- 
tions ■ A detailed description of MiniBooNE-specific 
modifications to NUANCE can be found in Ref. 
so only the most important details are reproduced here. 
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NUANCE uses a relativistic Fermi gas model to simu- 
late the carbon nucleus. The model is parametrized by 
a Fermi momentum of 220 ± 30 MeV/c and a binding 
energy of 34 ± 9 MeV, which are determined via electron 
scattering data [19j . 

Although the details of the CC7r + simulation are not 
essential to the extraction of the measured cross sec- 
tions, comparisons of the measurements with the de- 
fault MiniBooNE prediction will be shown in the re- 
sults section. In NUANCE, resonantly produced CC7r + 
events are simulated using the Rein-Sehgal model [2(| 
with M A = 1.10 ± 0.28 GeV/c 2 determined from ex- 
ternal experimental data [l]HL4 |. The model in NU- 
ANCE is further modified to include non-isotropic A- 
resonance decays according to Ref. [2(J. Pauli blocking 
is additionally accounted for in the decay of the reso- 
nance by requiring the momentum of the decay nucleon 
to be larger than the Fermi momentum. Coherently pro- 
duced CC7r + events are described using Rein-Sehgal [21] 
with M A = 1.03 ± 0.28 GeV/c 2 , and with the overall 
cross section rescaled by 0.65 to reproduce a prior Mini- 
BooNE measurement of coherent pion production in the 
NC channel [Hj]. 

The largest backgrounds to the CC7r + sample result 
from CCQE and CC multi-7r events. CCQE interac- 
tions are simulated using non-dipole vector form fac- 



tors [23j, a nonzero pseudoscalar form factor [24j, and 
a dipole axial-vector form factor with Ma — 1-23 ± 
0.08 GeV/c 2 [25| along with an additional Pauli blocking 
rcscaling, re = 1.02 ±0.02, as measured from MiniBooNE 
CCQE data [2||. These values for Ma and re are differ- 
ent from those recently reported in fl8j j but were chosen 
as they were extracted using the same, default model for 
resonance production in NUANCE as assumed here. An 
additional 10% normalization uncertainty is assumed to 
account for the differences between the relativistic Fermi 
gas and more modern nuclear models. Multipion pro- 
duction processes are modeled in NUANCE assuming 
M A = 1-30 ± 0.52 GeV/c 2 such that the sum of the ex- 
clusive CC channels reproduces CC inclusive data. 

For non-coherent scattering, NUANCE assumes neu- 
trino interactions take place on a single nucleon within 
the nucleus. The resulting particles (including reso- 
nances, nucleons, pions, etc.) can experience final-state 
interactions as they traverse the nuclear medium. For ex- 
ample, baryonic resonances can re-interact in the nucleus 
producing a pion-less final state with a probability of 20% 
for A++A and A°+A interactions, and 10% for A+++N 
and A - + N. An uncertainty of 100% is assumed for all 
four interaction probabilities 18]. Pions can also rescat- 
ter before exiting the target nucleus. Intranuclear pion 
absorption and charge-exchange processes are assigned 
uncertainties of 25% and 30%, respectively, based on ex- 
isting pion-carbon data |26l428j . Since the signal for this 
analysis is defined in terms of the particle content of the 
post-nuclear final state, the measurement uncertainty is 
not significantly affected by the uncertainties in intranu- 
clear pion rescattering. 



II. EVENT RECONSTRUCTION 

Neutrino events are reconstructed based on the charge 
and time recorded in each of the hit PMTs in the main 
tank volume. For a given set of seven initial track param- 
eters - energy, direction (9 and 4>), position (x, y, and z), 
and time - both a charge and a time probability distri- 
bution function (PDF) are produced for each hit PMT. 
The product of these PDFs evaluated at the measured 
charge and time values give a likelihood function, 

N nnhH N hit 

£(x) = J] P^unhit;x) n^(hit;x)/fe;x)/fe;x), 
t=i j=i 

(1) 

where Nhu(N un hit) is the number of hit(unhit) PMTs in 
the event and 7^ (hit; x)(7^(unhit; x)) is the probability 
that PMT i will be hit(unhit) for a track specified by x. 
The charge and time PDFs for x (/(g^jx) and f(tj\x), 
respectively) are evaluated at the measured charge, qj, 
and time, tj, in PMT j. The best set of parameters, 
x, are those that maximize the likelihood function. A 
complete description of the MiniBooNE extended-track 
reconstruction method is given in Ref. [29) . 



A. Pion Reconstruction 

To properly reconstruct CC7r + events, a mechanism for 
separating muons from charged pions is required. Unlike 
the case of separating muons and electrons, muons and 
pions propagate and emit Cherenkov radiation in a very 
similar manner due to their similar masses. The main 
differences are due to the hadronic interactions experi- 
enced by pions. When such a hadronic interaction takes 
place, the pion experiences an abrupt change in direction. 
Since the nuclear debris created in these interactions is 
generally well below Cherenkov threshold, the only de- 
tectable prompt light is produced by the "kinked" pion 
trajectory. 

To exploit these kinked pion trajectories, the straight- 
track fit hypothesis has been generalized to include four 
new parameters that characterize a kinked track. The 
length of the upstream portion of the track is determined 
by the energy lost prior to the kink point, AE up . The 
additional pion energy lost to the hadronic system during 
the interaction, AEkink , is also allowed to vary during the 
fit. Finally, the independent direction of the downstream 
track segment is characterized by two angles, Odown and 

4*down • 



1. Particle Identification Performance 

The ability of the kinked pion fitter to separate muons 
from pions is displayed in Fig. [TJ The peak in the pion 
likelihood ratio distribution is shifted relative to that of 
muons. More significantly, the pion distribution has a 
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FIG. 1: The straight muon and kinked pion likelihood ratios 
are shown for Monte Carlo muons (red) and pions (black) 
generated with full hadronic interactions and decays. The 
particles were generated from a flat kinetic energy distribution 
ranging from 50 to 450 MeV to approximate the true pion 
energy spectrum of CC7r + events. There is separation in the 
muon and pion peaks, and a large excess of pion events is seen 
along the high side tail. 



FIG. 2: The fractional pion kinetic energy reconstruction 
bias from the straight and kinked pion fits to Monte- Carlo- 
generated single pion events is shown. The low-energy shoul- 
der is significantly reduced in the kinked fitter, and rather 
than being 10% low, as is the case with the straight fitter, the 
peak from the kinked fitter is centered at zero. 



much larger tail of events that extend away from the 
muon portion of the likelihood ratio. These are events 
where kinked trajectories occurred and were successfully 
found by the fitter. 

The [ij~K separation provided by the kinked pion fit- 
ter is not as clean as the fi/e separation. There is no 
single value of the likelihood ratio at which a cut could 
be placed that would reject a large population of muons 
while retaining a significant fraction of pions. The goal 
of this analysis, however, is to reconstruct events with 
both a muon and a pion present, and to determine the 
identity of each track. In that case, the separation power 
indicated by Fig. Q] is doubled due to the presence of the 
second track. 



2. Kinematics Reconstruction Performance 

Although the main motivation for developing a kinked- 
track fitter was to provide a means for separating muons 
from charged pions, the improved modeling of pion tra- 
jectories results in superior event reconstruction as well. 
The fractional energy reconstruction bias (i.e. the ratio 
of the fit/true difference to the true value) from both 
the straight and kinked pion fitters is shown in Fig. [5] 
The straight pion fitter reconstructs pion energies 10% 
low, whereas the kinked fitter reconstruction bias peaks 



at zero. In addition, the "shoulder" just below the peak, 
where the reconstructed energy underestimates the true 
pion energy, is reduced by the kinked pion fitter. The 
two-dimensional plot of the fractional energy reconstruc- 
tion bias versus the true energy in Fig. [3] shows that 
the shoulder comes from higher energy pions that can 
produce multi-kink events and cause larger pion energy 
losses at each kink. The pion direction reconstruction 
is also significantly improved with the kinked fitter, as 
shown in Fig. [4] The event population in the first few 
bins of the angle between the reconstructed and true di- 
rections are nearly doubled in the kinked fitter relative 
to the straight fitter. 



B. CCtt Fit 

With the ability to reconstruct charged pions, a full 
CC7r + fitter is formed by simultaneously fitting for a 
straight muon and a kinked pion track. A CCn + fit has 
14 parameters: a common vertex (4 parameters), the ini- 
tial energy and direction of both the muon and pion (6 
parameters), and the additional kinked-track parameters 
for the pion (4 parameters). Just as in the kinked pion 
fitter, the predicted charges from all track segments (up- 
stream pion, downstream pion, and muon) are summed 
to get the total predicted charge for each PMT. 
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FIG. 3: The fractional pion kinetic energy reconstruction bias 
is plotted versus the true energy for Monte-Carlo-generated 
single pion events. The low-fit-energy shoulder is caused by 
higher energy pions. 



1. CCVr + Fit Seeding 

Each CC7r + event is assumed to have three Cherenkov 
rings from the three-track segments in the event: the up- 
stream pion segment, the downstream pion segment, and 
the muon track. Each ring is found in succession using 
intermediate two- and three-track likelihood functions. 
The three tracks are then pieced together in several dif- 
ferent pairings to create the kinked pion track and the 
straight muon track. The pairing that produces the best 
likelihood is used to seed the CC7r + fitter. 

The first of the three rings is found by performing a 
one-track fit. The results of the fit are frozen in place, and 
a scan for a second track is performed over 100 equally 
spaced directions. At each scan point, a two-track like- 
lihood function is evaluated, and the configuration that 
gives the best likelihood value is used to seed a full two- 
track fit. This process is then repeated by freezing the 
result of the two-track fit and scanning for a third track 
using a three-track likelihood function. The result of the 
scan seeds a full three-track fit. 

Once three tracks are found that characterize the 
Cherenkov rings in the event, they are combined to form 
a straight muon and a kinked pion. Only pairings where 
the downstream pion track has a lower energy than the 
upstream pion track are allowed. This reduces the num- 
ber of possible groupings to three. Each of these three 
seeds is passed to the full CC7T + fitter, and the fit that 
produces the best likelihood is chosen. 
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FIG. 4: The angle between the reconstructed and true pion 
directions is shown for both straight and kinked pion fits to 
Monte-Carlo-generated single pion events. The population in 
the lowest few bins where the properly reconstructed events 
lie is nearly twice as large for the kinked fitter. 



2. Fit Results 

The fractional kinetic energy reconstruction bias distri- 
butions for the muon and pion tracks are given in Fig.sJS] 
and[5J respectively. The muon kinetic energy has a small 
tail at low reconstructed energy due to fx/ir mispairing. 
The reconstructed pion kinetic energy has the same low- 
energy shoulder from high energy particles seen in the 
pion-only fit in Fig. [5] In addition, the fitter tends to 
place about 5% too much energy in the muon track at 
the expense of the pion track. 

While the pion energy fit is more accurate at low 
track energies, the opposite is true for the reconstructed 
pion direction. The track direction reconstruction re- 
lies on the existence of a well-defined Cherenkov ring 
from the upstream pion track segment. At MiniBooNE 
neutrino energies, 16% of the generated pion kinetic en- 
ergy spectrum lies below 70 MeV. This corresponds to 
an above-Cherenkov propagation distance of less than 
10 cm, which is often insufficient to determine the di- 
rection. For comparison, fewer than 1% of muons are 
generated below 70 MeV. In addition, the pion inelastic 
interaction length in mineral oil in the energy range of in- 
terest is approximately 1 m [26| . which means 10% of all 
pions will have upstream segments shorter than 10 cm. 

The ability of the fitter to correctly reconstruct both 
the muon and the pion directions is shown in Fig. [7] 
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FIG. 5: The fractional muon kinetic energy reconstruction 
bias is shown for all Monte Carlo signal events, and for cor- 
rectly paired signal events after all cuts other than the m„-+ m 
cut. Most of the low- fit-energy tail is due to events where the 
pion was misidentified as the muon. 



FIG. 6: The fractional pion kinetic energy reconstruction bias 
is shown for all Monte Carlo signal events, and for correctly 
paired signal events after all cuts other than the m^+N cut. 
The low-reconstructed-energy shoulder from the pion-only fit 
in Fig. [2] is seen here as well. 
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The reconstructed angle between the muon and pion is 
plotted against the larger of the two reconstructed/true 
angles. A perfect fitter would place all events in the 
lowest column where both reconstructed/true angles are 
zero. Fig. [7] shows that those bins contain the largest 
event population. The other significant event population 
is along the diagonal of the plot. These are events where 
the fitter has misidentified the muon as a pion and vice 
versa. In such cases, the angle between the true and 
reconstructed directions of both the muon and pion will 
be near the reconstructed fx/n angle. This is the first time 
charged pions have been tracked and their kinematics 
measured in a Chcrcnkov detector. 



3. Neutrino Energy 

With reconstructed energies and directions for both 
the muon and pion, the energy of the incident neutrino 
can be determined. Assuming the target nucleon is at 
rest and the remaining, unmeasured final-state parti- 
cle is a nucleon, the neutrino energy is specified by 4- 
momentum conservation, 

_ ml + ml - 2m N (E t _ t + E n ) + 2p^ ■ p w 

2 (Ef, + £L - |p M | cosfl^ - Ip^l cos6 Ui7r - m N ) ' 

(2) 




Max(6(Muon fit, true),6(Pion fit,true)) (Radians) 



FIG. 7: The reconstructed angle between the muon and pion 
directions for Monte Carlo CC7r + events is shown versus the 
larger (i.e. worse reconstructed) of the two reconstructed/true 
angles: 0(/x rec , Utrue) and 6(n rec ,ntrue)- The bins in the left- 
most columns represent events where both tracks have been 
properly reconstructed. Events in which the tracks have been 
misidentified appear along the diagonal. 
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FIG. 8: The neutrino energy reconstruction bias is plotted 
against the true neutrino energy for Monte Carlo generated 
CC7r + events. The reconstructed and true values are well 
correlated over the entire energy spectrum. 
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FIG. 9: The reconstructed n + N mass distribution is shown 
for both the data and the Monte Carlo simulation with full 
systematic uncertainties. The MC distribution has been nor- 
malized to the data. The signal and background components 
of the Monte Carlo distribution are also shown. At Mini- 
B00NE energies, the majority of CC7r + events come from de- 
cays of the A(1232) resonance. 
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where m„ E x , p x , and |p x | are the mass, energy, 4- 
momentum, and 3-momentum magnitude of particle x in 
the detector frame, and 9 uli (6 un ) is the angle between 
the directions of the muon(pion) and the neutrino. The 
neutrino direction is determined by the event vertex loca- 
tion and the mean neutrino emission point from the beam 
Monte Carlo prediction, although the large distance be- 
tween the beam and the detector means this angle is 
never larger than one degree. The comparison between 
reconstructed and true neutrino energy is given in Fig. [8] 
The resolution is 13.5% over most of the sensitive range, 
with a slight increase at the highest energies. 

The decreased pion angular resolution at lower pion 
energies has little impact on the neutrino energy recon- 
struction since the neutrino energy calculation becomes 
less sensitive to the reconstructed pion direction as the 
pion energy is reduced. In addition, events with misiden- 
tified tracks that are otherwise well-reconstructed will 
produce nearly the same neutrino energy, since muons 
and pions have similar masses. 



4- Invariant mass of the hadronic system 

By making the aforementioned assumptions required 
to calculate the neutrino energy, the kinematics of the 
interaction are fully specified. Previous attempts to mea- 
sure CC7r + interactions by reconstructing only the muon 
required the additional assumption that the recoiling par- 
ticle was an on-shell A baryon [3(|. Since the width of 
the A resonance is about 10% of its mass, this assump- 



tion results in an irreducible contribution to the neutrino 
energy resolution. By measuring the pion kinematics, the 
A mass constraint can be removed. 

The absence of a A mass constraint also means that the 
7r + + N invariant mass, which is dominated by the A res- 
onance, can be measured. Fig. |H] shows the reconstructed 
tt + + N mass, and a breakdown of the background com- 
position is given in Fig. [TO] The CCQE background fea- 
tures a sharp peak near threshold. CCQE interactions 
typically do not produce a pion, and the fitter correctly 
assigns very little kinetic energy to the hadronic system 
in these events. 



5. Momentum Transfer 

The final variable measured in this analysis is the 4- 
momentum transfer, q, from the leptonic current to the 
hadronic portion of the decay, which is characterized by 
its relativistic invariant, 

Q 2 = -{p^~p v ) 2 . (3) 

Since Q 2 is a property of the exchanged W boson, it is 
completely specified by the change in the leptonic cur- 
rent. However, this also means that, unlike the neutrino 
energy calculation, the reconstructed Q 2 distribution is 
quite sensitive to [i/ir misidentification. Fig. [TT] shows 
the fractional error in the reconstructed Q 2 distribution, 
normalized in columns of true Q 2 . Most of the columns 
peak near zero, but at high Q 2 , a second population of 
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FIG. 10: The reconstructed Monte Carlo ir + +N background 
distribution is divided into CCQE background events, and all 
other backgrounds. Since the CCQE events do not contain a 
pion, they are peaked near threshold (m, + Jnjv). 



FIG. 11: The Q 2 (fit-true) /true distribution for Monte Carlo 
events is plotted versus the true Q 2 . Each column of true Q 2 
has been normalized to unity. The reconstruction is able to 
determine the true Q 2 over the full range of the measurement. 
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events appears in which the fit underestimates the true 
Q 2 . These events are dominated by a high energy muon 
that has been misidentified as a pion, giving the impres- 
sion that most of the neutrino momentum was trans- 
ferred to the hadronic system. After all analysis cuts, 
the Q 2 resolution is 18% below 0.3 GeV 2 /c 4 and 20% 
below 1 GeV 2 /c 4 . 

III. EVENT SELECTION 

MiniBooNE events are recorded in a 19.2 us time win- 
dow beginning 4.6 fj,s prior to the arrival of the 1.6 /is 
beam pulse. The event time window is further subdi- 
vided into groups of PMT hits called subevents. The hit 
times in each subevent must have no more than two gaps 
longer than 10 ns and no gaps longer than 20 ns. 

The first subevent in each event consists of the final- 
state particles produced in the neutrino interaction. 
Muons and pions can stop in the detector and decay to 
produce Michel electrons. These Michel electrons result 
in additional subevents that provide a simple and pow- 
erful tool for separating neutrino event types. CCQE 
events, which contain a muon in the final state, most of- 
ten produce two subevents. Since CC7r + events are more 
likely to produce two Michel electrons from the decays 
of the final-state muon and pion, events are required to 
have three subevents. 

To remove backgrounds from cosmic rays, fewer than 
six hits are allowed in the veto region in all three 
subevents. The effect of this cut on the second subevent 



is shown in Fig. 1121 If a cosmic muon enters the tank 
before the beginning of the event time window, it can 
stop and decay within the event time window to simu- 
late a neutrino interaction. These events are removed by 
requiring a minimum number of PMT hits in the main 
tank in the first subevent. If this cut is placed at 200 
hits, more than 99.9% of beam unrelated backgrounds 
are removed fl6j . For the present analysis, the tank 
hits requirement has been relaxed to 175 hits since it 
is unusual for a Michel electron event to produce three 
subevents. The second and third subevents are required 
to have between 20 and 200 hits in the main tank to 
accept subevents from Michel electrons. 

To remove events that occur close to the edge of the de- 
tector, the muon and pion tracks are required to travel at 
least 150 cm before reaching the wall. Particle trajecto- 
ries that begin near the edge of the tank and are directed 
toward the tank wall are poorly reconstructed since they 
are detected by a small number of PMTs. Conversely, 
events that occur just inside the tank wall but consist of 
tracks that all point toward the interior of the tank are 
well reconstructed. An illustration of this cut is given in 
Fig.s[Hand[TU 

A final cut on the ir + N mass (m^iy) is used to elim- 
inated events where the final-state particles are misiden- 
tified. The fitter misreconstructs the muon as a pion, 
and vice versa, 21.4% of the time. In cases where a high 
energy muon is misreconstructed as a pion, the recon- 
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FIG. 12: The tank hits distribution in data is shown for the 
second subevent before and after requiring fewer than 6 hits in 
the veto. The Michel electron peak is mostly unaffected, while 
the large tail from entering comic rays is mostly removed. 



FIG. 13: The distance between the tank wall and event vertex 
along the muon trajectory is shown for both data and Monte 
Carlo events. All other cuts have been applied and the error 
bars include only data statistics. The cut on this distribution 
removes all events below 150 cm. 



FIG. 14: The distance between the tank wall and event vertex 
along the pion trajectory is shown for both data and Monte 
Carlo events. All other cuts have been applied and the error 
bars include only data statistics. The cut on this distribution 
removes all events below 150 cm. 
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structed kinematics give a spuriously high m^+N- The 
relationship between the generated and reconstructed 
771^+ jv is shown in Fig. [15] Beyond reconstructed masses 
of 1350 MeV/c 2 , the population of misreconstructed 
events begins to dominate, so a cut is implemented to 
remove these events. Fig. [TO] shows the improvement 



in the reconstructed muon kinetic energy resulting from 
the mn+N cut. The properly reconstructed events are 
mostly retained while the tail at low reconstructed en- 
ergy is greatly reduced. 

A well-matched event is defined in terms of the an- 
gles between true and reconstructed tracks. There are 
four such angles: f?(true /i, fit fi), f?(true 7r, fit tt), 
#(true [i, fit 7r), f?(true 7r, fit fi). If the minimum of 
these four angles is between a true track and its cor- 
responding reconstructed track, the event is said to be 
well-matched. This means that an event containing two 
properly reconstructed tracks is well-matched if the fit- 
ter correctly identifies the muon and pion tracks. Events 
where the direction of one of the tracks is misrecon- 
structed (e.g. one track is below Cherenkov threshold) 
are still considered to be well-matched if the angle be- 
tween the measured track and its corresponding true 
track is small. The fraction of well-matched events in- 
creases from 78.6% to 88.0% with the introduction of the 
m^+N < 1350 MeV/c 2 cut. 

After all cuts, 48,322 events are seen in the data with 
an overall signal efficiency of 12.7%, and a purity of 
90.0%. The background contributions are labeled accord- 
ing to the particles produced in the initial neutrino inter- 
action (i.e. prior to any final-state effects), rather than 
the final state emerging from the nucleus. The largest 
backgrounds (listed by percentage of the total sample) 
are from CC multi-7r events (3.1%), CCQE events (2.7%), 
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FIG. 15: The Monte Carlo m p i+jv distribution shows a cor- 
relation between the reconstructed and true distributions at 
low mass. At high reconstructed mass, the distribution is 
dominated by events with a high energy muon misidentificd 
as a pion. A cut is placed at 1350 MeV/c 2 to remove these 
events. 



and CC7r + events (1.3%) in which the pion content of the 
final state is altered via nuclear effects and are therefore 
not considered signal events. The complete lists of both 
signal and background compositions are given in Tables U 
and[TTl respectively. 



IV. CROSS SECTION MEASUREMENTS 

The observable CC-7r + analysis includes measurements 
of the cross section in terms of several kinematic vari- 
ables. The integrated cross section has been measured as 
a function of neutrino energy. The other one-dimensional 
measurements are differential cross sections as a function 
of the muon and pion energies and Q 2 . Since these one- 
dimensional measurements are necessarily averaged over 
the full neutrino energy spectrum, a corresponding two- 
dimensional measurement of each variable is made in bins 
of neutrino energy. In addition, the energy and direction 
are measured together for both the muon and pion in two 
double-differential cross section measurements. 



A. Data Unfolding 
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FIG. 16: The fractional muon kinetic energy reconstruction 
bias is shown for Monte Carlo CC7r + events before and after 
the m-K+N cut with all other cuts applied. The lower plot 
shows the fraction of events that pass the cut in each bin. 
The tail at low-fit-energy is significantly reduced. 



Because of biases and imperfect resolution in the event 
reconstruction, the event distributions measured in the 
data do not perfectly reflect the underlying true distri- 
butions. For example, as shown in Fig. [31 the pion energy 
reconstructed by the CC7r + fitter is systematically high 
at low true energy, and falls below the true energy as the 
energy increases. Since such biases are modeled in the 
Monte Carlo simulation, it is possible to "unfold" these 
bin-migration effects. 

The Monte Carlo bin-migration matrix for a given vari- 
able, v, is constructed by forming a two-dimensional his- 
togram of the reconstructed value of v versus the true 
value, and normalizing each true column to unity as illus- 
trated in Fig. 1171 Each element, Fji, in the bin-migration 
matrix represents the probability that an event generated 
with a value of v in bin i will be reconstructed in bin j. 
By definition, 



3 



(4) 



where ATj n * is the reconstructed distribution and Ti is the 
true distribution. 

To produce an unfolding matrix, the reconstructed- 
versus-true histogram is, instead, normalized in recon- 
structed rows. This produces a matrix, My, that per- 
forms the inverse operation of Eq. 2) This method, pro- 
posed by D'Agostini, avoids the problem of incorporating 
large statistical variance in the unfolding matrix, which 
is often an issue with inverting the bin-migration ma- 
trix [32j . The unfolding matrix can be extended to two- 
dimensional distributions in a straight forward way by 
arbitrarily ordering each of the two-dimensional bins and 
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TABLE I: The signal efficiency and purity are shown after each successive analysis cut. The efficiency is given relative to all 
interactions that occur within the detector, including the outer veto region; if restricted to events produced within 100 cm of 
the tank center, the effiency rises to 27%. A full list of efficiency and purity by NUANCE event type is given in Ref. [311 ]. 



Cut # Description Effic.(%) Purity (%) 



1 


no cuts 


100 


18.7 


2 


1st subevent, tank hits > 175 and veto hits < 6 


47.5 


23.6 


3 


number of subevents = 3 


23.2 


60.5 


4 


2nd and 3rd subevents: tank hits < 200 and veto hits < 6 


19.6 


86.3 


5 


muon and pion distance to wall > 150 cm 


16.9 


87.2 


6 


ty + N mass < 1350 MeV/c 2 


12.7 


90.0 



TABLE II: The composition of the background after all analysis cuts is given in terms of the NUANCE interaction mode 
(i.e. before final-state interactions) . The total background is 10% of the final event sample. The background acceptance and 
contamination after each cut is given in Ref. [311 ] . 



NUANCE Interaction 


Initial State Description 


Background Fraction (%) 


CC multi-7r 


1 muon, >1 pion, and 1 nucleon 


31.4 


CCQE 


1 muon and 1 proton 


26.8 


CC7T + 


1 muon, 1 7r + , and <1 nucleon 


13.4 


CCmesonB 


1 muon, 1 non-7r meson, and 1 baryon 


8.5 


CCDIS 


1 muon and multiple hadrons 


6.6 


CC^° 


1 muon, 1 7T°, and 1 proton 


5.7 


V 


all v interactions 


1.1 


other 




6.5 



repeating the same process used in the one-dimensional 
case. 

With the introduction of the unfolding matrix, My, 
the full expression for the differential cross section can 
be written as 

da Y^M^-B,) 

where vi is the variable to be measured in true bin i, Dj 
is the measured data distribution in reconstructed bin 
j, Bj is the predicted background distribution, q is the 
efficiency, Avi is the width of bin i, N tar g is the number 
of target molecules in the fiducial volume, and $ is the 
integrated flux in units of neutrinos per unit area, as 
described in Sec. IIVD1 

B. Unfolding Bias 

Although the use of My rather than F^ 1 avoids the 
statistical variance issues involved with matrix inversion, 
it does introduce some bias. In general, unfolding proce- 
dures often require the introduction of some amount of 
bias in order to reduce the statistical variances associated 
with matrix inversion such that the overall uncertainty 
is reduced 32]. Since the bin-migration matrix, Fy, is 
normalized in columns of the true distribution, Fij and 
F~j l are fully independent of the true Monte Carlo distri- 
bution, and are therefore unbiased transformations. The 
My matrix is normalized in reconstructed rows, which 
means any change to the shape of the true distribu- 
tion slice within a reconstructed bin (i.e. changes to the 



Bayesian prior probabilities as described in Ref. (32|) will 
result in the reconstructed events in that bin being as- 
signed to the true bins in different proportions. 

To quantify the size of the unfolding bias, an iterative 
technique is used. The background-subtracted, unfolded 
data provide an inferred true distribution as described in 
Section [IV Al Each Monte Carlo event is then assigned a 
weight given by the binned ratio of inferred true data to 
the true Monte Carlo simulation. Using these weights, a 
new reconstructed-versus-true histogram is created from 
which a new My unfolding matrix is formed. The data 
distribution is unfolded again using the new unfolding 
matrix and the processes is repeated. 

Successive iterations of the inferred data distribution 
have two distinct features in both the one- and two- 
dimensional cases. The first is that the largest excursion 
relative to the uniterated inferred distribution is in the 
first iteration. The other is that each successive iteration 
oscillates about an intermediate preferred value, which 
is a convolution of the true underlying distribution and 
any systematic biases in the unfolding matrix. The am- 
plitude of these oscillations decreases as the number of 
iterations increases. The size of the largest systematic 
variation, the first iteration, is taken as the systematic 
uncertainty. 



C. Efficiency Correction 

After the data are unfolded, the inferred true data dis- 
tribution is corrected for events lost due to data selection 
cuts and detector inefficiency. The numerator of the ef- 
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FIG. 17: The reconstructed- versus-true (top), bin- migration 
(middle), and unfolding (bottom) matrices are shown for 
muon kinetic energy. The bin-migration matrix is formed 
by normalizing the true columns of the Monte Carlo 
reconstructed-versus-true matrix to unity, while the unfold- 
ing matrix is formed by normalizing the reconstructed rows. 



ficiency correction is the true distribution of all Monte 
Carlo events that pass the cuts. The denominator is the 
generated Monte Carlo distribution, formed before any 
cuts are imposed. The ratio of these two distributions 
gives the fraction of events in a particular bin that sur- 
vive the analysis cuts, 



N, 



true after cuts 



^generate 



(6) 



The efficiency is insensitive to changes in the underlying 
physical parameters used to produce the generated dis- 
tribution. If any portion of the generated distribution is 
incorrectly enhanced, a proportional effect should be seen 
in the true distribution and thus cancel in the efficiency. 

Monte Carlo events are generated out to a radius of 
610.6 cm to include all neutrino interactions in the main 
tank, the veto region, and the tank wall. Since the mea- 
surement being performed is a neutrino cross section in 
mineral oil, all other materials must be excluded in form- 
ing the generated Monte Carlo distribution. To avoid the 
PMTs and, in particular, the material voids inside the 
PMTs, the efficiency denominator is formed from a sub- 
set of these events generated within a radius of 550 cm. 

Nearly all of the events generated outside of 575 cm 
are removed by the veto hits cut; however, there will be 
a population of events generated between 550 cm and 
575 cm that pass all cuts, particularly in the upstream 
portion of the tank. The extra contribution from events 
in the 550-575 cm shell are in good agreement in the 
reconstructed data and Monte Carlo vertex distributions. 

The number of interaction targets in the cross section 
formula, N tar g, corresponds to the definition of the gen- 
erated volume used in the efficiency calculation. To ex- 
tract the number of targets from the volume, the only 
experiment-dependent quantity needed is the oil density, 
which has been measured to be 0.845 ±0.001 g/cm 3 [la ]. 
Since the cross section only depends on the relative 
amount of each atomic species, the interaction target is 
chosen to be an average single unit on the hydrocarbon 
chain, CH2.08- 



D. Flux Factor 

The flux factor, in Eq. [5] takes the form of either a 
distribution in neutrino energy or a single value, depend- 
ing on the type of cross section measurement being per- 
formed. For measurements binned in E Vl the flux factor 
is the number of incident neutrinos per unit area in each 
measured bin. These measurements are flux-averaged 
over the width of each bin. In the differential and double- 
differential cross section measurements, the flux factor is 
the fully integrated neutrino flux, and the cross section is 
flux-averaged over the entire neutrino energy spectrum. 
The total integrated flux factor for the MiniBooNE 
experiment, normalized to protons on target (POT), is 
5.19 x 10~ 19 (^/cm 2 /POT). 

Flux-averaged differential cross section measurements 
implicitly contain the shape of the neutrino energy spec- 
trum, which must also be reported to fully specify the 
results. To mitigate the dependence on the experiment- 
dependent shape of the energy spectrum, each of the dif- 
ferential cross section measurements has also been per- 
formed in bins of neutrino energy. 
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FIG. 18: The reconstructed muon kinetic energy is plotted 
in each of the 67 optical model multisims. The central value 
Monte Carlo distribution (red) is overlayed for comparison. 



E. Systematic Uncertainties 



The systematic uncertainties are grouped by error 
sources. Each source is a set of correlated uncertainties 
from a particular stage of the simulation. The parameters 
of each error source are related by an error matrix from 
which a set of correlated parameter variations, called a 
"multisim," can be drawn. Each new set of parameters 
produces a systematically varied version of any recon- 
structed distribution. The spread in the reconstructed 
multisim distributions is used to calculate the total sys- 
tematic uncertainty. 

There are two distinct types of systematic variations. 
Some systematics, such as the flux and cross section 
uncertainties, only affect the probability with which an 
event will occur. For this type of uncertainty, a systemat- 
ically varied distribution can be produced by reweighting 
the central value Monte Carlo distribution. Each event 
is multiplied by the ratio of the event probability calcu- 
lated with the systematically varied set of parameters to 
the central value event probability. 



1. Error Matrices 

The uncertainties in the measured cross sections are 
described by an error matrix that characterizes the cor- 
related uncertainties in the measured values in each bin. 
For each error source, an error matrix, Ef° urce , is calcu- 
lated from the bin population differences in the multisims 
compared to the central value, 

1 N 

E source _ — (n m ,i - n C vs) (n m ,j - ncvj) , (7) 

m— 1 

where N is the number of multisims, n m j is the number 
of events in bin i of multisim m, and ncv,i is the num- 
ber of events in bin i in the central value Monte Carlo 
simulation. Once an error matrix has been calculated for 
each source, the total error matrix is given by summing 
each component matrix. 

Since the statistical fluctuations in a particular bin are 
unrelated to the fluctuations in any other bin, statisti- 
cal error matrices are always diagonal. By design, these 
uncertainties are built into the optical model error ma- 
trix since each optical model multisim was constructed 
to have the same number of events as the data. Un- 
fortunately, this also has the effect of adding statistical 
fluctuations to the off-diagonal terms. As more optical 
model multisims are incorporated into the calculation of 
the error matrix, the size of these spurious fluctuations 
is decreased. The fluctuations are also smaller for bins 
with significant event populations. For this reason, cross 
section results will only be reported for bins with at least 
25 unfolded data events. In the one-dimensional cross 
section measurements, the event populations are large 
enough that this effect is negligible. The reweighting 
multisims do not suffer from this effect. 

To evaluate the systematic uncertainties in the cross 
section, the full cross section calculation procedure out- 
lined in Eq. [S] is performed in each multisim. The multi- 
sim distributions replace the corresponding central value 
distributions in the calculation. The reconstructed data 
distribution remains the same, but multisim distributions 
are used for the unfolding matrix, the background pre- 
diction, and the signal efficiency. The formula for the 
differential cross section with multisim dependent quan- 
tities explicitly identified is 



The other type of systematic variation affects the mea- 
sured properties of the event after it is produced, such 
as the amount of light generated as a function of wave- 
length and the propagation of the light through the oil. 
In general, these variations cannot be accomplished via 
reweighting. Instead, these errors are determined using 
67 data-sized Monte Carlo simulations that are gener- 
ated using parameter draws from the optical model error 
matrix. A plot of these optical model multisims in muon 
kinetic energy is shown in Fig. 1181 



dv 



(8) 



where m is the multisim index. 

The results presented in the remainder of this report 
will list only the diagonal errors on each bin. The full 
error matrices for each distribution are quite large, and 
in the case of the two-dimensional measurements, they 
can contain over one million elements. 
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FIG. 19: The energy spectrum is shown (top) along with 
the fractional flux errors (bottom). The fractional errors are 
cumulative such that each successive error band includes a 
quadrature sum of the previous errors such that the outer 
band gives the total error. A numerical description of the 
energy spectrum is given in Table IXIXI 



tion, and the resulting fit parameters are systematically 
varied according to the error matrix returned by the fit. 
The covariance of these spline function variations with 
respect to the Sanford-Wang central value are used to 
form the beam tt + production error matrix as described 
in Ref. [18j . This produces an uncertainty given by the 
HARP data errors in regions where the Sanford-Wang 
parametrization agrees with the spline function, and in- 
creases the errors in regions where they differ. 

The portion of the flux relevant to observable CC7r + 
interactions occurs at neutrino energies larger than about 
400 MeV. For the peak neutrino energies (0.5-1 GeV), 
the beam ir + flux uncertainties are around 8%. At higher 
neutrino energies, the beam ir + uncertainty grows to as 
large as 25%. The beam ir + fluctuations also exhibit 
very distinct features in shape. The residual effect of the 
wiggling behavior of the spline fit to the HARP data is 
apparent. In particular, the low-energy flux exhibits very 
large systematic excursions since there are no HARP data 
to constrain the fits in that region. 

The other two flux related error sources have a smaller 
effect on the total uncertainty. The "Beam" uncertain- 
ties, which contain all systematic effects involving the 
proton beam and horn, are generally around the 5% level 
below 1 GeV and then expand at higher neutrino ener- 
gies. The main contributor at high energies is the horn 
current skin depth uncertainty, which causes a ~15% ef- 
fect between 1.5 and 2.5 GeV. The K + production un- 
certainties are mostly irrelevant for this analysis. K + 
mesons become the dominant source of v ^ production at 
2.3 GeV, and the uncertainties become dominant at neu- 
trino energies greater than 2.5 GeV, where the flux is 
very small. 

For the measurements not involving neutrino energy, 
the cross section calculation is only affected by the uncer- 
tainty in the integrated flux. The size of these variations 
for each of the flux error sources is given in Table Mil 



2. Flux Uncertainties 

The predicted neutrino spectrum is modified to cor- 
respond to the systematically varied flux parameters in 
each multisim. The diagonal uncertainties of the flux 
variations are shown in Fig. [121 

The largest of the flux errors is the uncertainty in 
the beam tt + production in proton-beryllium interac- 
tions. The differential cross section for these interactions 
is parametrized using an empirical model from Sanford 
and Wang [33]. The parameters are determined in a si- 
multaneous fit to data from the HARP and E910 exper- 
iments [13, H|[ . The flux prediction is not tuned in any 
way to data from MiniBooNE. As described in Ref. (la ], 
the shape of the Sanford-Wang parametrization is not 
fully compatible with the data, but it is still used in the 
Monte Carlo simulation to provide physical constraints 
such as driving the cross section to zero at vanishing pion 
momentum. To determine the uncertainties in the pion 
production, the HARP data are fit with a spline func- 



TABLE III: The uncertainties in the integrated flux are given 
for each of the flux error sources. 



Error Source 


Integrated Flux Uncertainty 


n + 


10.4% 


Beam 


4.1% 


K+ 


0.4% 



Results 



The observable v^GG^ 



cross section on a CH2.08 tar- 
get has been measured in a variety of forms: a total cross 
section as a function of neutrino energy (Fig. |20|) . a dif- 
ferential cross section in Q 2 (Fig. |2"TT) . differential cross 
sections in the kinetic energy of the muon (Fig. [22]) and 
the pion (Fig. |2"3"]) . double-differential cross sections in the 
muon kinetic energy and angle (Fig. I24[) and the pion ki- 
netic energy and angle (Fig. 1251) . and two-dimensional 
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measurements of each of the differential cross sections in 
bins of neutrino energy to provide results independent 
of the MiniBooNE energy spectrum (Fig.s [26l [27l and 
|2"5)) . This is the first time model- independent differential 
cross sections have been provided for the muon and pion 
kinematics in these interactions. 

The binning for each of the one-dimensional distribu- 
tions has been chosen such that the true Monte Carlo 
prediction in each bin exceeds 250 events after all cuts. 
The one-dimensional bin sizes are used for the two- 
dimensional measurements as well to retain sufficient pre- 
cision in the most interesting regions of phase space. This 
results in several bins with very small numbers of pre- 
dicted events. The data-sized optical model multisims 
produce unreliable uncertainties in bins with small event 
populations, and therefore results will only be reported 
for bins that contain at least 25 inferred true data events. 

TABLE IV: The uncertainties in the total, flux-averaged cross 
section are given for the dominant error sources. 



Error Source 


Cross Section Uncertainty 


Beam ty + Production 


9.2% 


v Cross Sections 


8.2% 


Proton Beam and Horn 


4.3% 


Optical Model 


1.5% 


Other 


<3% 



The uncertainties from the most significant error 
sources in the total cross section, averaged over the neu- 
trino energy spectrum, are shown in Table IIVI In each 
of the one-dimensional differential cross section measure- 
ments, the two largest sources of uncertainty are the 
beam 7r + production and the neutrino interaction cross 
sections. The 7r + production uncertainties in the flux- 
averaged results are dominated by the large uncertainties 
at low neutrino energy. Since the low-energy region has 
relatively little impact on the measurements binned in 
neutrino energy, the beam ir + production uncertainties 
are is significantly lower and generally remain below 10% 
except at the highest neutrino energies. 

The largest effects in the cross section uncertainties are 
pion absorption and charge- exchange interactions that 
take place after the pion has left the target nucleus. If 
the pion is absorbed, it will not produce a Michel electron 
and the event will fail the three-subevent requirement; 
therefore, pion absorption and charge- exchange interac- 
tions will directly affect the cut efficiencies. A 50% un- 
certainty is assigned to the pion charge-exchange cross 
section and a 35% uncertainty is assigned to pion ab- 
sorption based on the agreement between the GCALOR 
Monte Carlo simulation [36| and external data [26l - |28| . 
The remainder of the cross section uncertainty is due to 
variations in the interaction cross sections of each back- 
ground process. 

Since the cross section measurements in Q 2 and neu- 
trino energy include contributions from the incident neu- 
trino, they must be unfolded back to the initial neutrino 
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FIG. 20: The <j{E v ) measurement is shown with cumulative 
systematic errors. The absolutely normalized Monte Carlo 
prediction is shown for comparison. The bottom plot shows 
the fractional uncertainties and the ratio of the Monte Carlo 
prediction to the measurement. 



interaction, and are therefore dependent on the modeling 
of nuclear effects. In particular, additional uncertainties 
in the kinematics of the target nucleons are absorbed in 
these results. Conversely, the measurements in the muon 
and pion kinematic variables are properties of only the 
final, post-nuclear state, and are therefore largely insen- 
sitive to nuclear model uncertainties. 

Finally, for most of the measured phase space, the un- 
certainty due to unfolding bias is negligible; however, it 
becomes significant at low Q 2 in both the one- and two- 
dimensional measurements. This particular region has 
two features that generally make unfolding difficult. The 
first is that the shape is rapidly changing, which strongly 
affects bin migration. Also, this is a region where the 
shapes in the data and Monte Carlo simulation signifi- 
cantly disagree, which increases the probability that the 
shape of the true Q 2 distributions within each recon- 
structed bin are incorrect. Despite these features, the 
unfolding uncertainty is still not the dominant system- 
atic effect, and is of comparable size only in the few bins 
at low Q 2 which are susceptible to these effects. 



16 



Error Bands 

• MiniBooNE Measurement 
| Total Uncertainty 
MC Prediction 




3. 40 



Error Bands 

• MiniBooNE Measurement 
Total Uncertainty 
MC Prediction 



1000 1200 1400 
Q 2 (MeV 2 /c 4 ) 




800 1000 1200 1400 
Muon Kinetic Energy (MeV) 



Error Bands 

• MiniBooNE Measurement 
Total Uncertainty 
MC/Measurement 




1.5 
1.4 
1.3 
1.2 
1.1 
1 

0.9 
0.8 
0.7 
0.6 
0.5, 



Error Bands 

• MiniBooNE Measurement 
Total Uncertainty 
MC/Measurement 




800 1000 1200 1400 
Muon Kinetic Energy (MeV) 



FIG. 21: The da/d(Q 2 ) measurement is shown with cumu- 
lative systematic errors. The absolutely normalized Monte 
Carlo prediction is shown for comparison. The bottom plot 
shows the fractional uncertainties and the ratio of the Monte 
Carlo prediction to the measurement. 



FIG. 22: The da/d(K E,j) measurement is shown with cumu- 
lative systematic errors. The absolutely normalized Monte 
Carlo prediction is shown for comparison. The bottom plot 
shows the fractional uncertainties and the ratio of the Monte 
Carlo prediction to the measurement. 



V. CONCLUSION 

Results have been presented for the observable CC7r + 
cross section as a function of several fundamental kine- 
matic variables. Of these results, the cross sections mea- 
sured as a function of the neutrino energy are the least 
experiment dependent, since the predicted neutrino flux 
has been accounted for separately in each bin. Previous 
measurements of the neutrino energy cross section below 
2 GeV have uncertainties larger than 20%, upon which 
the present results provide a significant improvement. 

The present measurement is on average 23% higher 
than the NUANCE prediction for the observable cross 
section, which is compatible with some previous results, 
most notably Ref. [14| . However, a direct comparison to 
past data is difficult since previous cross section measure- 
ments in this energy range were conducted on hydrogen 
or deuterium. Since MiniBooNE employs a nuclear tar- 
get, it is unclear if the source of discrepancy seen in the 
current results lies in the single nucleon cross section, or 
whether it is due to nuclear effects. 



The remaining kinematic CC7r + cross sections have not 
been reported previously. Since the single- and double- 
differential cross section measurements are necessarily 
averaged over the shape of the neutrino flux prediction 
given in Fig. 1191 each measurement of a final-state kine- 
matic quantity has also been measured in bins of neutrino 
energy to remove the dependence on the MiniBooNE en- 
ergy spectrum. The integrated one-dimensional measure- 
ments have also been included due to the familiarity of 
many in the community with flux-averaged results. This 
is the most complete set of information that has ever 
been available for CC7r + on nuclear targets. 
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FIG. 23: The dcr /d(KEn) measurement is shown with cumu- 
lative systematic errors. The absolutely normalized Monte 
Carlo prediction is shown for comparison. The bottom plot 
shows the fractional uncertainties and the ratio of the Monte 
Carlo prediction to the measurement. 
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FIG. 24: The measured d 2 a/d(cos(e i _ l> „))d(KE l _ l ) values are 
shown (top) along with the total fractional uncertainties (mid- 
dle). Empty bins indicate regions where no measurement has 
been made. The Monte Carlo predicted cross section is shown 
for comparison (bottom). 
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FIG. 25: The measured d 2 a /d(cos(6 7T , v ))d(K E^) values are 
shown (top) along with the total fractional uncertainties (mid- 
dle). Empty bins indicate regions where no measurement has 
been made. The Monte Carlo predicted cross section is shown 
for comparison (bottom). 



FIG. 26: The measured da(E v ) / d(Q 2 ) values are shown (top) 
along with the total fractional uncertainties (middle) . Empty 
bins indicate regions where no measurement has been made. 
The Monte Carlo predicted cross section is shown for com- 
parison (bottom). 
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FIG. 27: The measured da(E v ) / 'd(KE^) values are shown 
(top) along with the total fractional uncertainties (middle). 
Empty bins indicate regions where no measurement has been 
made. The Monte Carlo predicted cross section is shown for 
comparison (bottom). 



FIG. 28: The measured do(E v )/d{KE 1T ) values are shown 
(top) along with the total fractional uncertainties (middle). 
Empty bins indicate regions where no measurement has been 
made. The Monte Carlo predicted cross section is shown for 
comparison (bottom). 
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Appendix A: Cross Section Tables 

Tables IVllXVIIII give the numerical values of each mea- 
sured cross section and Table IXlXl gives the MiniBooNE 
flux prediction. 

TABLE V: The o{E v ) results from Fig.E0]are given with the 
total uncertainty. Each row is a bin of neutrino energy (MeV) 
labeled according to its low edge. 

Bin (MeV) Result (lCT^cm^) 
500 6.1 ± 0.8 



600 


11.5 


± 


1.3 


650 


15.8 


± 


1.7 


700 


19.6 


± 


2.1 


750 


24.1 


± 


2.5 


800 


28.3 


± 


2.9 


850 


32.5 


± 


3.3 


900 


37.3 


± 


3.9 


950 


41.6 


± 


4.4 


1000 


46.6 


± 


5.0 


1050 


49.7 


± 


5.6 


1100 


52.9 


± 


6.1 


1150 


56.3 


± 


7.0 


1200 


59.1 


± 


7.6 


1250 


62.3 


± 


8.4 


1300 


66.7 


± 


9.4 



1350 70.3 ± 10.5 

1400 72.3 ± 11.4 

1450 77.6 ± 12.8 

1500 80.8 ± 13.9 

1550 83.7 ± 14.9 

1600 86.4 ± 15.9 

1650 88.5 ± 16.9 

1700 93.0 ± 18.5 

1750 92.6 ± 19.3 

1800 97.1 ± 21.4 

1900 99.2 ± 23.7 



TABLE VI: The dafd(KE^) results from Fig. [23] are given 
with the total uncertainty. Each row is a bin of pion kinetic 
energy (MeV) labeled according to its low edge. 



Bin (MeV) Result (10~ 4i cm7MeV) 



o 


3.6 ± 0.5 


50 


8.9 ± 1.2 


75 


10.8 ± 1.4 


100 


10.8 ± 1.2 


125 


9.7 ± 1.2 


150 


8.6 ± 1.1 


175 


7.6 ± 1.1 


200 


7.0 ± 1.0 


225 


6.1 ± 1.0 


250 


5.3 ± 0.9 


275 


4.5 ± 0.8 


300 


3.9 ± 0.7 


325 


3.2 ± 0.6 


350 


2.7 ± 0.5 


375 


2.3 ± 0.4 



TABLE VII: The da/d{Q 2 ) results from Fig. [2T] are given 
with the total uncertainty. Each row is a bin of Q 2 (GeV 2 /c 4 ) 
labeled according to its low edge. 

"Bin (GeV7c 4 ) Result (10~ 45 cm^c 4 /MeV 7 )" 



0.00 


36.2 ± 6.7 


0.05 


55.8 ± 7.8 


0.10 


56.1 ± 7.5 


0.15 


53.4 ± 7.0 


0.20 


47.8 ± 6.4 


0.25 


43.2 ± 5.9 


0.30 


38.8 ± 5.3 


0.35 


34.1 ± 4.8 


0.40 


30.3 ± 4.3 


0.45 


26.1 ± 3.7 


0.50 


22.9 ± 3.3 


0.55 


19.8 ± 2.9 


0.60 


17.3 ± 2.5 


0.65 


14.8 ± 2.1 


0.70 


12.9 ± 1.9 


0.75 


11.0 ± 1.7 


0.80 


9.3 ± 1.4 


0.85 


8.0 ± 1.3 


0.90 


6.8 ± 1.1 


0.95 


5.3 ± 1.0 


1.05 


3.9 ± 0.7 


1.15 


2.5 ± 0.5 


1.30 


1.4 ± 0.4 
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TABLE VIII: The da/d(KE,,) results from Fig.[22]are given 
with the total uncertainty. Each row is a bin of muon kinetic 
energy (MeV) labeled according to its low edge. 

"Bin (MeV) Result (l(T^cm7MeV) 



o 


23 1 +32 


50 


34.4 ± 5.0 


100 


3Q + 5 8 


150 


41.6 ± 6.2 


200 


41.2 ± 6.1 


250 


3Q + ^ 7 


300 


37.3 ± 5.4 


350 


3"i + 4 9 


400 


32.0 ± 4.4 


450 


28 8 + 3 8 


500 


2fi 1 +35 


550 


23 9 + 3 1 


DUU 


90 £ -1- 9 & 


650 


18.6 ± 2.5 


700 


15.9 ± 2.1 


750 


14.2 ± 2.0 


800 


12.3 ± 1.8 


850 


10.7 ± 1.6 


900 


9.3 ± 1.4 


950 


7.7 ± 1.3 


1000 


5.9 ± 1.0 


1100 


4.3 ± 0.8 


1200 


2.2 ± 0.6 
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TABLE IX: The da(E v )/d(Q 2 ) results from Fig. [26] are shown (10" 44 cm 2 c 4 /MeV 2 ). Each bin is labeled according to its low 
edge. The columns are bins of neutrino energy (MeV) and the rows are bins of Q 2 (GeV 2 /c 4 ). Empty bins indicate regions 
where no measurement has been made. 
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1.4 


2.4 
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3.3 
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5.4 


5.2 


5.0 


5.9 
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0.55 
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2.4 


3.0 


3.8 


4.1 


4.6 


4.7 


5.2 


5.1 


5.7 


5.6 


6.1 


5.8 
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1.8 


2.4 
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3.6 
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4.8 


5.4 


5.6 


5.8 


5.8 


6.7 


6.6 


6.6 
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6.2 


7.5 


7.3 


7.1 


0.45 










1.5 


2.3 


3.1 


4.0 


4.7 


5.2 


5.7 


5.5 


6.1 


6.5 


6.4 


6.6 


6.9 


7.7 


7.1 


7.4 


6.6 


7.9 


8.1 


0.40 








1.3 


2.5 


3.2 


4.2 


4.8 


5.1 


6.1 


6.5 


6.7 


6.7 


7.1 


7.1 


7.3 


7.5 


7.1 


8.1 


7.9 


8.6 


8.4 


9.3 


0.35 






1.1 


2.0 


3.1 


4.0 


4.8 


5.7 


6.4 


6.8 


6.8 


7.0 


7.4 


7.8 


7.6 


8.2 


8.1 


8.1 


8.7 


9.1 


8.7 


8.9 


8.8 


0.30 






1.9 


3.1 


4.1 


5.1 


6.0 


6.1 


7.4 


7.4 


7.4 


8.0 


8.0 


8.1 


8.6 


9.1 


9.2 


9.0 


8.7 


8.7 


8.9 


8.7 


9.8 


0.25 




1.8 


3.0 


4.1 


5.3 


5.9 


6.7 


6.8 


7.5 


7.7 


8.2 


8.2 


8.1 


8.5 


9.0 


9.0 


9.8 


9.4 


9.4 


8.8 


9.5 


9.7 


10.1 


0.20 


1.0 


2.9 


4.2 


5.2 


6.0 


6.7 


6.9 


7.6 


8.0 


8.8 


8.5 


8.8 


8.7 


9.3 


9.2 


9.5 


9.6 


9.5 


9.9 


10.3 


10.2 


11.5 


11.2 


0.15 


1.9 


4.1 


5.4 


6.3 


7.1 


7.4 


7.5 


8.3 


8.7 


9.1 


9.0 


9.3 


9.2 


9.6 


9.4 


9.6 


9.6 


9.9 


10.6 


11.6 


10.8 


11.2 


10.3 


0.10 


3.1 


5.2 


6.6 


7.0 


7.4 


7.6 


7.9 


8.7 


8.2 


9.2 


9.3 


9.4 


9.5 


9.1 


9.0 


9.5 


9.2 


9.5 


10.4 


9.4 


10.5 


10.2 


9.8 


0.05 


4.3 


6.0 


6.5 


7.0 


7.3 


7.4 


7.6 


8.0 


7.9 


8.0 


8.0 


8.3 


8.6 


8.2 


8.6 


8.2 


8.4 


7.8 


8.6 


9.3 


9.3 


9.6 


10.0 


0.00 


1.9 


2.7 


3.3 


3.6 


3.9 


4.3 


4.6 


4.8 


5.2 


5.4 


5.4 


5.2 


5.8 


5.9 


6.0 


6.4 


6.6 


6.1 


6.3 


6.7 


7.2 







TABLE X: The percent uncertainty of the da(E v )/d(Q 2 ) results from Fig. [26] is shown. Each bin is labeled according to its 
low edge. The columns are bins of neutrino energy (MeV) and the rows are bins of Q 2 (GeV 2 /c 4 ). Empty bins indicate regions 
where no measurement has been made. 



Bin 


500 


600 


650 


700 


750 


800 


850 


900 


950 


1000 


1050 


1100 


1150 


1200 


1250 


1300 


1350 


1400 


1450 


1500 


1550 


1600 


1650 


1.15 


































25.7 


26.5 


34.4 










1.05 




























27.9 


33.9 


26.0 


21.6 


22.1 


33.5 










0.95 


























22.7 


23.5 


20.8 


21.9 


21.5 


21.5 


20.6 


25.1 


22.0 


21.9 




0.90 


























21.8 


25.5 


22.6 


19.8 


21.7 


21.1 


22.8 










0.85 






















24.1 


20.0 


22.1 


18.8 


20.0 


22.1 


19.6 


21.3 


20.6 










0.80 




















26.4 


22.7 


20.2 


18.9 


17.7 


18.5 


19.5 


17.9 


21.8 


21.8 










0.75 




















22.2 


22.5 


17.3 


16.1 


16.8 


16.4 


18.9 


18.9 


20.5 


20.0 


23.0 








0.70 


















21.0 


15.9 


19.4 


17.2 


15.7 


15.1 


15.6 


18.3 


18.7 


22.6 


20.9 


25.0 


23.7 






0.65 
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25.6 
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13.6 
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19.9 
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0.55 
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17.4 
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19.1 
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0.45 










19.7 
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13.5 


13.9 
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14.4 
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15.3 


16.4 
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27.1 


29.2 


0.40 
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15.9 


13.9 
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0.35 
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14.8 
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0.30 






13.0 
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14.5 
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0.25 
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14.8 


13.3 


11.2 
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11.1 
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12.3 
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12.7 
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14.3 
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0.15 


10.9 


10.6 


11.2 


11.0 


10.9 


11.1 
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12.0 
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11.9 
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13.7 


14.0 
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19.4 
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10.3 


11.2 


10.8 


10.9 
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11.9 


11.5 


12.1 


11.9 


12.0 


12.1 


13.6 


14.5 


14.8 


15.7 


16.1 


16.7 


17.6 


18.4 


19.0 


21.5 


23.1 


0.05 


12.2 


12.8 


12.2 


13.9 


13.8 


13.3 


12.8 


11.8 


11.9 


11.7 


12.8 


13.0 


13.7 


14.8 


15.8 


17.5 


18.4 


21.1 


20.4 


19.1 


21.4 


20.9 


23.9 


0.00 


20.0 


18.6 


19.1 


18.9 


17.7 


17.2 


16.2 


14.9 


14.8 


13.6 


15.8 


15.3 


14.9 


14.7 


15.3 


17.3 


20.9 


20.7 


21.9 


24.0 


21.2 
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TABLE XI: The da{E„)/d(KE^) results from Fig. [27] are shown (10" 41 cm 2 /MeV). Each bin is labeled according to its low 
edge. The columns are bins of neutrino energy (MeV) and the rows are bins of muon kinetic energy (MeV). Empty bins indicate 
regions where no measurement has been made. 



Bin 


500 


600 


650 


700 


750 


800 


850 


900 


950 


1000 


1050 


1100 


1150 


1200 


1250 


1300 


1350 


1400 


1450 


1500 


1550 


1600 


1650 


1200 














































1.9 


1100 










































4.1 


6.3 


8.0 


1000 




































1.8 


3.6 


5.9 


8.0 


8.9 


9.0 


950 


































2.8 


4.6 


6.9 


8.6 


9.4 


9.5 


9.6 


900 
































2.5 


4.9 


6.9 


8.8 


9.4 


9.8 


9.1 


9.5 


850 






























2.5 


4.2 


6.8 


8.4 


9.1 


10.3 


9.7 


9.5 


9.1 


800 


























1.2 


2.4 


4.7 


7.1 


7.9 


9.5 


9.4 


9.6 


9.5 


8.5 


9.2 


750 
























1.3 


2.7 


4.8 


7.0 


8.5 


9.5 


9.4 


9.5 


8.3 


8.3 


9.1 


9.5 


700 






















1.1 


2.4 


4.8 


7.0 


8.1 


8.9 


8.8 


9.1 


9.1 


8.5 


8.8 


8.6 


7.5 


650 




















1.2 


2.5 


4.6 


6.9 


8.4 


9.7 


9.7 


9.9 


9.2 


8.9 


8.0 


8.1 


7.3 


7.4 


600 


















1.2 


2.5 


4.8 


7.0 


8.7 


9.1 


9.0 


9.1 


9.2 


9.3 


8.4 


8.4 


7.8 


7.9 




550 
















1.2 


2.4 


4.7 


6.7 


8.3 


9.1 


9.4 


8.9 


9.8 


8.8 


7.6 


8.9 


7.9 


8.1 


7.1 




500 














1.2 


2.6 


4.5 


7.1 


8.4 


9.4 


9.2 


9.3 


9.1 


8.9 


8.9 


8.0 


8.2 


8.0 


7.1 


7.4 




450 












1.1 


2.4 


4.5 


7.1 


8.3 


9.1 


9.8 


9.7 


9.0 


8.7 


8.6 


8.2 


8.9 


7.4 


7.7 


7.6 


6.3 




400 










1.0 


2.2 


4.3 


7.0 


8.7 


9.6 


9.9 


9.4 


8.9 


9.1 


9.1 


8.4 


7.9 


7.5 


7.9 


7.0 


6.5 


6.2 




350 








1.0 


2.1 


4.3 


6.8 


8.3 


9.8 


10.4 


9.9 


9.5 


9.1 


8.4 


8.9 


8.4 


7.6 


7.1 


7.0 


7.5 


6.1 






300 




0.3 


1.0 


2.1 


4.3 


6.2 


7.9 


9.3 


9.5 


10.2 


10.2 


9.1 


9.4 


8.4 


8.4 


8.3 


7.3 


7.2 


7.3 


6.7 


6.7 






250 


0.2 


0.9 


2.0 


3.9 


6.2 


7.8 


8.5 


9.5 


9.8 


9.4 


9.4 


9.0 


8.5 


8.3 


7.6 


7.7 


6.7 


6.8 


6.6 


6.0 








200 


0.5 


1.8 


4.0 


6.0 


7.7 


8.5 


9.3 


9.5 


9.6 


9.3 


8.9 


8.3 


8.0 


7.5 


6.8 


6.4 


6.8 


6.0 


5.7 


5.8 








150 


1.3 


3.7 


5.7 


7.6 


8.5 


9.3 


9.3 


8.9 


8.5 


8.2 


7.5 


7.4 


6.7 


6.7 


6.2 


5.7 


5.4 


4.7 


4.9 










100 


2.2 


5.4 


7.1 


8.0 


8.5 


8.5 


7.9 


7.5 


6.5 


6.6 


6.1 


5.6 


5.4 


4.8 


4.5 


4.0 


4.1 














50 


4.0 


6.7 


7.5 


7.2 


7.0 


6.2 


5.6 


5.1 


4.5 


4.2 


3.9 


3.6 


3.4 


3.2 























4.5 


4.8 


4.6 


3.8 


3.5 


3.1 


2.7 


2.6 


2.4 


2.4 


2.2 



























TABLE XII: The percent uncertainty of the da(E^)/d(KEfj,) results from Fig. [27] is shown. Each bin is labeled according to 
its low edge. The columns are bins of neutrino energy (MeV) and the rows are bins of muon kinetic energy (MeV) . Empty bins 
indicate regions where no measurement has been made. 



Bin 


500 


600 


650 


700 


750 


800 


850 


900 


950 


1000 


1050 


1100 


1150 


1200 


1250 


1300 


1350 


1400 


1450 


1500 


1550 


1600 


1650 


1200 














































28.6 


1100 










































22.9 


22.9 


24.8 


1000 




































21.9 


22.0 


19.8 


21.2 


23.3 


23.8 


950 


































22.2 


19.6 


19.7 


20.3 


20.9 


22.1 


26.8 


900 
































19.2 


19.5 


20.1 


19.5 


20.1 


22.4 


23.3 


25.2 


850 






























20.7 


19.0 


17.9 


17.9 


18.8 


19.8 


23.3 


23.7 


25.9 


800 


























21.2 


20.6 


19.9 


16.8 


17.4 


17.7 


19.2 


21.5 


22.0 


24.5 


25.8 


750 
























18.3 


16.4 


15.1 


15.6 


15.9 


16.1 


18.1 


19.1 


21.1 


20.5 


26.1 


30.5 


700 






















20.4 


18.0 


15.6 


15.3 


15.8 


16.0 


16.5 


17.9 


18.5 


20.4 


21.3 


22.3 


28.9 


650 




















20.0 


18.4 


15.6 


15.1 


14.3 


15.3 


16.1 


18.0 


17.9 


19.6 


21.0 


21.9 


22.0 


25.8 


600 


















15.8 


14.9 


13.6 


12.9 


13.9 


15.6 


14.9 


15.5 


17.2 


19.2 


20.1 


20.8 


22.6 


25.6 




550 
















15.5 


14.4 


14.1 


13.5 


13.7 


13.7 


15.1 


15.2 


16.4 


17.0 


19.1 


20.9 


21.3 


24.7 


22.9 




500 














15.0 


12.3 


12.5 


12.7 


11.9 


13.4 


13.8 


14.4 


15.6 


17.0 


17.8 


18.2 


20.8 


25.2 


24.9 


26.4 




450 












17.0 


13.7 


12.3 


11.8 


10.8 


11.7 


12.3 


13.5 


14.2 


14.8 


15.5 


16.4 


19.2 


19.2 


21.7 


23.0 


22.8 




400 










18.8 


16.3 


13.7 


11.8 


11.5 


11.5 


11.8 


12.3 


13.8 


13.7 


16.5 


16.8 


17.0 


18.7 


20.9 


20.5 


23.5 


23.9 




350 








20.3 


17.4 


14.5 


12.0 


10.9 


11.0 


12.0 


12.5 


13.2 


13.8 


14.8 


15.8 


16.8 


18.8 


19.6 


18.9 


20.4 


21.6 






300 




21.5 


16.9 


15.4 


14.0 


12.3 


11.4 


10.9 


11.1 


12.3 


13.9 


15.3 


15.4 


15.5 


16.1 


20.0 


19.1 


21.2 


21.3 


21.3 


24.1 






250 


22.5 


16.5 


13.6 


12.6 


12.0 


11.9 


10.4 


11.0 


11.7 


12.4 


13.3 


14.7 


15.6 


16.3 


18.4 


17.4 


19.4 


22.9 


22.8 


23.0 








200 


18.2 


14.1 


11.7 


10.4 


10.7 


10.8 


10.9 


12.3 


12.2 


12.7 


14.4 


15.0 


16.2 


17.7 


17.9 


20.0 


19.3 


19.4 


21.4 


24.2 








150 


14.6 


12.6 


10.8 


10.7 


11.0 


11.5 


11.3 


12.0 


13.4 


13.7 


15.4 


14.8 


15.7 


17.0 


18.9 


18.7 


20.3 


21.7 


20.6 










100 


14.2 


11.6 


10.8 


10.7 


10.6 


11.4 


11.9 


12.1 


13.0 


13.7 


14.3 


15.2 


16.4 


18.1 


17.7 


19.9 


20.5 














50 


13.3 


12.9 


11.2 


12.0 


11.9 


11.7 


13.8 


14.9 


13.9 


15.4 


16.7 


17.5 


17.9 


18.5 























12.8 


11.9 


13.5 


12.7 


14.3 


14.6 


15.6 


16.2 


16.9 


18.2 


18.4 
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TABLE XIII: The da(E„)/d(KE^) results from Fig. [28] are shown (10~ 41 cm 2 /MeV). Each bin is labeled according to its low 
edge. The columns are bins of neutrino energy (MeV) and the rows are bins of pion kinetic energy (MeV). Empty bins indicate 
regions where no measurement has been made. 



Bin 


500 


600 


650 


700 


750 


800 


850 


900 


950 


1000 


1050 


1100 


1150 


1200 


1250 


1300 


1350 


1400 


1450 


1500 


1550 


1600 


1650 


375 












1.6 


1.8 


2.3 


2.7 


3.5 


4.3 


4.8 


5.1 


5.2 


6.0 


6.8 


7.0 


7.5 


8.0 










350 










1.3 


1.8 


2.2 


3.0 


3.5 


4.4 


5.2 


5.4 


6.0 


6.9 


6.7 


7.2 


8.2 


8.0 


8.9 


9.5 








325 








1.1 


1.9 


2.6 


3.2 


4.0 


5.0 


5.2 


6.0 


7.0 


6.8 


7.5 


7.9 


8.6 


9.2 


10.1 


10.9 


10.9 


11.5 






300 






1.1 


1.6 


2.6 


3.7 


4.1 


4.8 


5.7 


6.6 


7.2 


8.0 


8.2 


8.7 


9.9 


10.9 


10.8 


11.0 


11.5 


12.3 


13.7 


12.8 




275 




0.9 


1.6 


2.7 


3.6 


4.4 


5.2 


6.4 


7.3 


7.8 


8.4 


8.7 


9.4 


10.4 


10.6 


11.7 


10.9 


12.0 


14.4 


13.5 


14.2 


15.8 


13.6 


250 




1.4 


2.4 


3.4 


4.6 


5.4 


6.3 


7.1 


8.2 


9.5 


10.0 


11.1 


11.6 


11.3 


11.8 


12.8 


12.8 


13.5 


14.3 


14.8 


16.5 


16.0 


16.7 


225 


0.7 


2.4 


3.3 


4.4 


5.5 


6.4 


7.5 


8.3 


9.5 


11.0 


11.2 


12.1 


12.8 


12.8 


13.4 


14.3 


15.4 


15.7 


15.3 


16.2 


17.9 


17.5 


17.2 


200 


1.4 


3.1 


4.2 


5.6 


6.9 


8.3 


9.2 


10.0 


11.3 


12.3 


12.0 


13.1 


14.9 


14.2 


15.0 


15.8 


16.5 


16.8 


17.8 


18.0 


17.7 


20.3 


18.8 


175 


1.9 


3.8 


5.2 


6.3 


8.1 


9.5 


10.0 


11.3 


12.3 


13.2 


13.1 


14.9 


15.5 


15.6 


16.0 


16.9 


18.2 


17.2 


19.9 


20.6 


19.1 


19.7 


22.0 


150 


2.6 


5.2 


6.9 


7.9 


8.8 


9.9 


11.4 


12.8 


13.6 


15.2 


16.3 


15.8 


16.5 


16.1 


17.7 


18.2 


18.8 


20.2 


18.9 


20.4 


22.1 


22.4 


19.8 


125 


3.4 


5.7 


7.2 


8.6 


10.5 


11.2 


12.4 


14.7 


15.9 


17.5 


17.6 


18.0 


18.4 


19.5 


19.7 


21.2 


20.6 


22.0 


21.2 


22.2 


24.0 


24.1 


24.4 


100 


3.8 


6.5 


8.9 


10.6 


11.7 


12.7 


14.5 


15.8 


17.1 


18.5 


19.0 


20.3 


20.2 


21.1 


21.9 


22.1 


22.6 


22.7 


25.0 


25.3 


25.3 


26.5 


26.7 


75 


4.2 


7.0 


9.0 


10.5 


12.3 


13.7 


14.9 


16.2 


16.8 


18.3 


19.7 


19.3 


18.9 


20.7 


20.8 


22.0 


22.9 


21.7 


23.3 


25.0 


24.1 


21.6 


24.8 


50 


3.7 


5.9 


7.4 


8.6 


9.9 


11.2 


12.4 


13.7 


14.1 


15.4 


14.7 


15.5 


16.3 


16.8 


18.0 


17.0 


17.8 


17.8 


18.8 


19.0 


18.9 


19.2 











2.9 


3.4 


3.8 


4.5 


5.1 


5.4 


5.8 


6.1 


6.6 


6.6 


7.0 


6.8 


7.0 


7.7 


8.1 


7.7 













TABLE XIV: The percent uncertainty of the da(E v )/d(KE 7r ) results from Fig. [28]is shown. Each bin is labeled according to 
its low edge. The columns are bins of neutrino energy (MeV) and the rows are bins of pion kinetic energy (MeV). Empty bins 
indicate regions where no measurement has been made. 



Bin 


500 


600 


650 


700 


750 


800 


850 


900 


950 


1000 


1050 


1100 


1150 


1200 


1250 


1300 


1350 


1400 


1450 


1500 


1550 


1600 


1650 


375 












22.5 


20.1 


20.9 


19.5 


20.8 


20.5 


20.8 


20.8 


20.6 


21.2 


20.3 


23.3 


23.2 


24.2 










350 










23.2 


20.2 


18.9 


19.3 


18.2 


19.0 


20.3 


19.3 


20.7 


19.9 


19.2 


21.1 


23.3 


23.0 


26.8 


25.8 








325 








20.6 


19.4 


18.8 


19.3 


18.8 


18.2 


18.4 


17.2 


18.6 


19.3 


19.9 


20.1 


20.2 


22.3 


21.3 


22.5 


24.3 


23.5 






300 






22.6 


19.0 


17.0 


18.2 


17.9 


17.6 


16.6 


16.8 


17.7 


17.6 


19.7 


20.0 


19.0 


19.9 


20.3 


21.9 


27.2 


21.9 


22.8 


23.2 




275 




23.0 


19.9 


17.6 


17.3 


17.3 


16.4 


17.2 


16.0 


17.0 


17.1 


17.5 


18.4 


18.4 


17.9 


18.9 


19.9 


19.7 


19.6 


24.3 


21.6 


23.0 


25.5 


250 




19.8 


18.1 


17.1 


15.6 


15.5 


15.1 


16.2 


16.0 


16.0 


16.2 


15.8 


17.5 


16.7 


17.9 


18.2 


18.7 


18.7 


21.7 


21.3 


24.6 


22.7 


26.3 


225 


19.4 


17.3 


16.7 


16.0 


15.3 


15.1 


15.0 


14.7 


14.9 


14.6 


14.3 


14.4 


16.7 


15.9 


18.4 


16.7 


17.2 


19.4 


22.8 


24.8 


19.4 


22.5 


24.9 


200 


17.0 


16.3 


15.6 


15.3 


13.7 


14.1 


14.1 


13.5 


13.7 


13.0 


13.5 


13.4 


15.7 


15.2 


15.9 


16.2 


18.5 


18.4 


18.8 


20.0 


21.2 


21.9 


23.4 


175 


16.6 


17.3 


15.9 


14.1 


12.9 


12.8 


12.5 


12.2 


12.2 


12.7 


13.3 


13.5 


13.9 


13.8 


14.8 


16.8 


16.2 


17.4 


18.6 


19.4 


19.9 


22.5 


24.4 


150 


15.6 


14.9 


13.5 


12.8 


12.5 


12.1 


11.5 


11.6 


11.9 


11.3 


12.5 


12.4 


13.3 


14.9 


14.7 


14.8 


16.2 


18.1 


19.5 


20.2 


20.6 


22.6 


21.3 


125 


16.5 


14.6 


11.6 


10.8 


11.5 


11.1 


10.5 


10.5 


10.2 


10.5 


11.4 


12.1 


13.1 


13.0 


14.2 


14.7 


16.1 


17.3 


17.7 


17.4 


19.9 


20.0 


24.4 


100 


14.4 


11.9 


10.9 


9.8 


10.0 


10.1 


9.2 


9.8 


9.7 


11.5 


10.8 


11.1 


11.9 


12.9 


13.2 


14.3 


15.4 


15.7 


17.2 


17.4 


18.3 


20.1 


21.2 


75 


12.0 


12.8 


10.4 


11.3 


11.6 


13.0 


12.0 


12.8 


12.7 


12.7 


14.0 


13.7 


13.2 


13.8 


14.4 


14.2 


15.3 


16.6 


15.9 


18.7 


18.6 


18.8 


21.5 


50 


14.5 


13.9 


11.9 


13.9 


14.3 


13.9 


15.4 


14.3 


13.3 


14.1 


14.4 


14.7 


16.7 


14.5 


13.9 


15.0 


17.1 


16.6 


15.9 


18.0 


18.5 


18.3 











15.3 


16.0 


15.9 


15.9 


16.1 


16.5 


12.9 


14.0 


14.3 


14.0 


14.9 


13.6 


13.8 


14.7 


14.2 


15.5 













25 



TABLE XV: The d 2 a/d{KE li )d{cos(9^)) results from Fig.[24]are shown (10" 42 cm 2 /MeV). Each bin is labeled according to 
its low edge. The columns are bins of muon kinetic energy (MeV) and the rows are bins of cos(muon, neutrino angle). Empty 
bins indicate regions where no measurement has been made. 



0111 


10U 




a a c 

u.yt> 


OO A 

16.0 


20.1 


a aa 
0.9U 


OA C 

30.0 


o c o 
30.8 


a o c 
0.O0 


Ol 7 

31. ( 


|)7 A 

3 1 .0 


A OA 

O.oU 


Ol *7 

31. ( 


OA 1 

39.1 


A 7C 
0. 10 


33. 1 


/I A A 

4U.0 


0.70 


37.8 


40.5 


0.65 


35.9 


41.7 


0.60 


33.5 


43.3 


0.55 


35.0 


41.7 


0.50 


32.8 


38.3 


0.45 


32.8 


37.8 


0.40 


32.5 


38.1 


0.35 


36.9 


40.5 


0.30 


31.1 


34.7 


0.25 


30.2 


31.8 


0.20 


27.2 


31.6 


0.15 


26.2 


26.8 


0.10 


25.3 


25.3 


0.05 


24.7 


25.4 


0.00 


24.5 


23.1 


-0.05 


24.6 


18.6 


-0.10 


21.3 


17.9 


-0.20 


21.0 


16.7 


-0.30 


17.5 


12.3 


-0.40 


14.7 


10.2 


-0.50 


13.1 


8.3 


-0.60 


11.6 


5.7 


-0.80 


8.5 


3.0 


-1.00 


3.9 





250 300 350 400 450 500 550 600 650 700 750 800 850 900 950 1000 1100 1200 



32.1 
41.3 
44.8 
43.4 



46.8 
47.0 
46.1 



29.9 
30.2 
27.1 
23.1 
21.3 
20.6 
17.6 
13.1 
12.4 
10.2 
5.7 
6.0 



38.5 


42.8 


47.9 


55.8 


58.2 


59.8 


63.7 


65.8 


69.3 


46.7 


55.0 


59.5 


64.3 


70.8 


70.3 


72.6 


73.3 


65.0 


50.9 


58.7 


70.2 


67.6 


67.5 


70.4 


71.4 


69.5 


56.0 


53.8 


63.6 


62.5 


64.7 


67.3 


61.7 


57.2 


55.3 


45.4 


53.3 


57.3 


58.6 


58.5 


60.4 


53.0 


44.9 


37.0 


33.1 


58.5 


57.3 


58.2 


53.8 


52.4 


42.4 


35.5 


28.8 


22.0 


50.6 


56.1 


49.2 


44.8 


38.4 


29.6 


23.5 


16.9 


14.1 


47.1 


51.6 


51.1 


40.6 


31.0 


25.0 


16.7 


11.5 


8.5 


49.4 


42.3 


37.3 


33.2 


24.1 


17.1 


11.1 






44.2 


39.5 


33.6 


26.1 


15.9 


11.8 








39.1 


32.8 


25.8 


17.9 


13.3 


8.8 








37.7 


28.9 


22.4 


14.9 


9.2 










31.6 


27.1 


17.0 


12.1 












28.8 


20.9 


12.3 














22.9 


17.9 


11.2 














22.3 


14.2 


10.6 














18.4 


9.4 
















15.6 


















13.3 


















10.3 



















71.5 65.9 58.5 

65.7 59.4 58.0 

52.7 49.3 42.3 
37.1 33.6 25.1 

25.8 18.7 13.7 
15.4 12.3 
10.0 



59.3 56.0 

50.8 46.4 
36.1 27.3 

19.9 15.5 



48.4 
38.3 
19.8 
9.0 



43.3 
28.1 
12.9 



27.7 
12.6 
5.7 



5.4 
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TABLE XVI: The percent uncertainty of the d 2 a/d(KE ti )d(cos(0 lilU )) results from Fig. [24] is shown. Each bin is labeled 
according to its low edge. The columns are bins of muon kinetic energy (MeV) and the rows are bins of cos(muon, neutrino 
angle). Empty bins indicate regions where no measurement has been made. 



Bin 


150 


200 


250 


300 


350 


400 


450 


500 


550 


600 


650 


700 


750 


0.95 


27.0 


25.9 


20.7 


19.0 


17.5 


17.1 


14.7 


16.2 


16.3 


14.5 


15.0 


13.9 


13.9 


0.90 


18.9 


18.1 


17.6 


16.7 


16.0 


14.6 


13.9 


15.0 


14.3 


14.5 


14.3 


15.4 


14.6 


0.85 


20.3 


16.9 


16.0 


17.0 


16.4 


16.0 


13.8 


14.7 


13.6 


14.8 


14.5 


13.8 


14.3 


0.80 


18.5 


16.6 


16.7 


15.1 


15.8 


15.9 


14.2 


14.5 


13.2 


14.0 


14.3 


14.7 


15.1 


0.75 


17.0 


16.3 


15.9 


14.9 


14.0 


14.5 


14.5 


14.1 


14.2 


15.6 


16.4 


15.9 


16.5 


0.70 


15.9 


15.6 


15.6 


15.8 


16.5 


15.9 


14.7 


15.0 


14.8 


15.0 


15.7 


20.2 


20.5 


0.65 


16.4 


15.4 


14.5 


15.7 


15.4 


14.8 


15.4 


15.0 


16.2 


16.5 


16.1 


20.1 


26.5 


0.60 


15.9 


15.1 


15.3 


16.1 


16.0 


16.2 


15.1 


15.4 


16.1 


20.7 


21.0 


23.9 




0.55 


16.4 


16.8 


16.4 


15.1 


14.7 


15.1 


17.2 


18.9 


17.8 


21.8 








0.50 


16.4 


17.4 


15.1 


17.1 


16.4 


17.1 


16.8 


19.1 


25.5 










0.45 


15.5 


17.8 


17.0 


16.3 


15.6 


17.8 


19.0 


17.6 


26.0 










0.40 


16.4 


16.8 


17.1 


17.0 


16.0 


17.2 


18.7 


26.3 












0.35 


19.1 


20.3 


17.7 


16.2 


17.3 


22.5 


25.2 














0.30 


18.1 


20.6 


17.2 


15.5 


18.1 


20.8 
















0.25 


17.1 


16.9 


16.2 


16.8 


18.5 


24.5 
















0.20 


17.1 


17.8 


18.0 


20.2 


19.2 


35.3 
















0.15 


17.9 


16.6 


16.9 


21.1 


22.2 


















0.10 


16.4 


17.1 


17.7 


21.1 




















0.05 


18.7 


18.5 


21.0 


24.0 




















0.00 


18.4 


21.1 


22.5 


27.8 




















-0.05 


20.1 


20.0 


20.8 






















-0.10 


21.8 


17.7 


25.6 






















-0.20 


19.1 


18.9 


26.1 


36.8 




















-0.30 


17.7 


21.6 


28.3 






















-0.40 


18.7 


20.4 


35.2 






















-0.50 


21.7 


27.5 
























-0.60 


24.0 


32.9 
























-0.80 


24.2 


29.5 
























-1.00 


25.4 



























800 850 900 950 1000 1100 1200 



15.2 16.3 
13.8 14.4 
15.4 17.6 
15.7 18.2 

17.3 21.8 
21.0 



15.9 15.8 

14.0 16.8 

18.1 19.4 
22.4 27.1 



15.8 
16.8 
20.3 
33.1 



15.6 
20.6 
24.6 



19.7 
25.0 
36.6 
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TABLE XVII: The d 2 a/d{KE 7r )d{cos{6^, lJ )) results from Fig. [25]are shown (10" 41 cm 2 /MeV). Each bin is labeled according 
to its low edge. The columns are bins of pion kinetic energy (MeV) and the rows are bins of cos(pion, neutrino angle). Empty 
bins indicate regions where no measurement has been made. 



Bin 


150 


175 


200 


225 


250 


275 


300 


325 


350 


375 


0.95 


7.9 


8.9 


10.3 


10.7 


10.0 


10.3 


9.4 


8.0 


6.8 


6.0 


0.90 


8.5 


9.6 


10.1 


9.9 


9.9 


9.1 


8.8 


7 5 


6.3 


5.4 


0.85 


9.7 


10.0 


10.8 


10.0 


9.2 


8.5 


7.5 


6.7 


5.5 


4.7 


0.80 


9.2 


9.8 


9.6 


9.2 


8.4 


7.5 


6.9 


5.9 


5.3 


4.4 


U. t o 


9.6 


9.0 


9.1 


8.6 


7.8 


7.2 


6.3 


5.3 


4.2 


4.1 


0.70 


9.2 


9.7 


8.9 


8.2 


7.4 


6.3 


5.4 


4.3 


3.9 


3.4 


n fi<i 


9.4 


8.7 


8.2 


7.6 


6.5 


5.8 


4.9 


4.1 


3.2 


2.9 


0.60 


9.3 


8.4 


8.4 


7.0 


5.9 


5.0 


3.9 


3.4 


3.0 


2.5 


0.55 


9.0 


7.7 


7.2 


6.1 


5.5 


4.6 


3.5 


3.1 


2.4 


2.2 


U.<JU 


8 4 

0.4: 


7 r i 


fi A 


"> fi 


u . u 


o.o 


'i 9 
o . z 


9 fi 


2.1 






7 ^ 


fi 7 


■7. O 


4 


4.2 


o . o 


9 7 
z . / 








n An 


7 ^ 


fi 9 

U. Z 


"l 4 


4.4 


O . - J 


9 fi 
z .u 


2.2 








0.35 


7.0 


5.8 


4.9 


3.6 


3.1 


2.5 










0.30 


6.8 


5.5 


4.3 


3.3 


2.5 












0.25 


6.1 


5.1 


3.6 


2.6 


2.3 












n on 
u.zu 


^ 7 
O. / 


1 9 




9 7 
Z. t 














0.15 


5.1 


3.7 


2.8 


1.9 














0.10 


4.5 


3.5 


2.3 


1.7 














0.05 


4.3 


2.9 


2.4 
















0.00 


3.7 


2.4 


2.0 
















-0.05 


3.1 


2.1 


















-0.10 


2.8 


1.8 


















-0.15 


2.3 




















-0.20 


2.1 




















-0.25 


2.0 





















TABLE XVIII: The percent uncertainty of the d 2 a / d(KE 7I )d{cos{8 7T ^)) results from Fig. [25] is shown. Each bin is labeled 
according to its low edge. The columns are bins of pion kinetic energy (MeV) and the rows are bins of cos(pion, neutrino 
angle). Empty bins indicate regions where no measurement has been made. 



Bin 


150 


175 


200 


225 


250 


275 


300 


325 


350 


375 


0.95 


16.5 


16.7 


17.3 


17.8 


19.7 


20.8 


21.3 


22.6 


22.1 


22.8 


0.90 


14.4 


16.1 


15.9 


17.4 


18.1 


19.2 


19.8 


20.1 


22.0 


21.4 


0.85 


12.8 


13.5 


14.8 


16.0 


17.4 


18.1 


19.0 


18.3 


19.1 


20.0 


0.80 


13.8 


14.3 


14.8 


15.6 


17.5 


19.8 


19.2 


19.5 


20.1 


19.8 


0.75 


13.5 


15.3 


15.9 


15.8 


17.3 


17.7 


19.6 


18.7 


21.5 


20.3 


0.70 


14.3 


16.0 


15.3 


17.1 


17.6 


17.8 


19.9 


18.0 


19.0 


19.8 


0.65 


14.0 


16.0 


15.5 


16.8 


16.8 


17.7 


17.1 


22.4 


19.2 


21.8 


0.60 


14.0 


16.9 


16.4 


17.6 


19.6 


19.7 


22.2 


20.3 


21.7 


23.3 


0.55 


15.7 


17.4 


20.4 


19.4 


18.6 


19.7 


19.8 


21.2 


32.3 


21.1 


0.50 


16.2 


17.9 


18.8 


18.1 


19.2 


19.0 


22.6 


23.5 


20.0 




0.45 


17.0 


18.2 


19.3 


19.4 


19.4 


23.1 


23.0 








0.40 


17.1 


16.7 


19.5 


23.2 


22.5 


20.8 


25.2 








0.35 


16.8 


16.3 


18.4 


21.7 


21.9 


23.8 










0.30 


17.8 


19.0 


19.7 


20.8 


22.2 












0.25 


16.9 


18.7 


21.3 


22.9 


22.0 












0.20 


18.3 


18.8 


21.9 


20.8 














0.15 


18.5 


19.7 


18.5 


20.9 














0.10 


17.3 


16.6 


17.8 


20.5 














0.05 


17.1 


16.4 


17.5 
















0.00 


17.4 


20.4 


18.4 
















-0.05 


16.8 


17.6 


















-0.10 


18.4 


20.1 


















-0.15 


21.5 




















-0.20 


18.9 




















-0.25 


20.4 
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E v bin 


U/j, flux 


E v bin 


flux 


E v bin 


flux 


(GeV) 


(^ M /POT/GeV/cm 2 ) 


(GeV) 


(^/POT/GeV/cm 2 ) 


(GeV) 


(^/POT/GeV/cm 2 ) 


0.00-0.05 


4.54 


X 


10" 11 


1.00-1.05 


3.35 


X 


io- 1 " 


2.00-2.05 


1.92 


X 


IO" 11 


0.05-0.10 


1.71 


X 


io- 10 


1.05-1.10 


3.12 


X 


io- 10 


2.05-2.10 


1.63 


X 


io- 11 


0.10-0.15 


2.22 


X 


io- 10 


1.10-1.15 


2.88 


X 


io- 10 


2.10-2.15 


1.39 


X 


io- 11 


0.15-0.20 


2.67 


X 


io- 10 


1.15-1.20 


2.64 


X 


10 -io 


2.15-2.20 


1.19 


X 


io- 11 


0.20-0.25 


3.32 


X 


io- 10 


1.20-1.25 


2.39 


X 


io- 10 


2.20-2.25 


1.03 


X 


io- 11 


0.25-0.30 


3.64 


X 


io- 10 


1.25-1.30 


2.14 


X 


io- 10 


2.25-2.30 


8.96 


X 


io- 12 


0.30-0.35 


3.89 


X 


io- 10 


1.30-1.35 


1.90 


X 


io- 10 


2.30-2.35 


7.87 


X 


io- 12 


0.35-0.40 


4.09 


X 


IO" 10 


1.35-1.40 


1.67 


X 


io- 10 


2.35-2.40 


7.00 


X 


io- 12 


0.40-0.45 


4.32 


X 


io- 10 


1.40-1.45 


1.46 


X 


1Q -io 


2.40-2.45 


6.30 


X 


io- 12 


0.45-0.50 


4.48 


X 


IO" 10 


1.45-1.50 


1.26 


X 


10 -io 


2.45-2.50 


5.73 


X 


io- 12 


0.50-0.55 


4.56 


X 


io- 10 


1.50-1.55 


1.08 


X 


io- 10 


2.50-2.55 


5.23 


X 


io- 12 


0.55-0.60 


4.58 


X 


io- 10 


1.55-1.60 


9.20 


X 


io- 11 


2.55-2.60 


4.82 


X 


io- 12 


0.60-0.65 


4.55 


X 


IO" 10 


1.60-1.65 


7.80 


X 


io- 11 


2.60-2.65 


4.55 


X 


io- 12 


0.65-0.70 


4.51 


X 


IO" 10 


1.65-1.70 


6.57 


X 


io- 11 


2.65-2.70 


4.22 


X 


io- 12 


0.70-0.75 


4.43 


X 


IO" 10 


1.70-1.75 


5.52 


X 


io- 11 


2.70-2.75 


3.99 


X 


io- 12 


0.75-0.80 


4.31 


X 


io- 10 


1.75-1.80 


4.62 


X 


io- 11 


2.75-2.80 


3.84 


X 


io- 12 


0.80-0.85 


4.16 


X 


IO" 10 


1.80-1.85 


3.86 


X 


io- 11 


2.80-2.85 


3.63 


X 


io- 12 


0.85-0.90 


3.98 


X 


IO" 10 


1.85-1.90 


3.23 


X 


io- 11 


2.85-2.90 


3.45 


X 


io- 12 


0.90-0.95 


3.79 


X 


io- 10 


1.90-1.95 


2.71 


X 


io- 11 


2.90-2.95 


3.33 


X 


io- 12 


0.95-1.00 


3.58 


X 


10 -io 


1.95-2.00 


2.28 


X 


10" 11 


2.95-3.00 


3.20 


X 


IO" 12 



TABLE XIX: Predicted i/ M flux at the MiniBooNE detector. 
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